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Using master equations we present an analytical solution of the time evolution of an entangled 
electron spin pair which can occupy 36 different quantum states in a double quantum dot nanostruc- 

oo ; 

^S) , ture. This solution is exact given a few realistic assumptions and takes into account relaxation and 

decoherence rates of the electron spins as phenomenological parameters. Our systematic method 
' of solving a large set of coupled differential equations is straightforward and can be used to obtain 

I ' analytical predictions of the quantum evolution of a large class of complex quantum systems, for 

, which until now commonly numerical solutions have been sought. 
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. I. INTRODUCTION 

o 

^ ■ Master equations are used to describe the quantum evolution of a physical system interacting with some "reservoir''^, 

and have been applied to a wide variety of physical systems, ranging from two-level atoms in the presence of light 
> • 

psj ' fields^ to solid-state nanostructurcs such as quantum dots and Joscphson junction devices'^. For simple systems, such 
^ ■ 

\^ as a two-level atom damped by a reservoir consisting of simple harmonic oscillators or an electron in a single or 
00 ■ 

. double quantum dot coupled to external leads", the set of master equations that describes the quantum dynamics 
of the system is small and its solution can be obtained analytically in a straightforward way. If the system is more 



involved, however, due to the presence of quite a few atomic levels or because the nanostructure is composed of 
Q various coherent parts, its quantum state space consists of a large number of quantum states with various coherent 

O • equations does not spring to the eye. Hence often a numerical solution is sought®. Understanding the quantum 
^ \ evolution of such "complex" quantum systems - where "complex" refers to a system which is described by a large 
number of coupled quantum states - has recently become increasingly important, in particular in fundamental research 
aimed at investigating the dynamic behavior of qubits, the basic building blocks for quantum computation^. A large 
theoretical and experimental effort in various fields, e.g. quantum optics, atomic physics and condensed matter physics, 
is presently directed towards investigating possibilities to use two-level systems such as polarized photons, cold atoms, 
electron spins and superconducting circuits as qubits and finding ways to couple these qubits together. In the latter 
three systems, one of the major questions involved is how the desired coherent evolution of the system will be affected 
by coupling to the environment, which is necessary to manipulate and measure the states of the qubits but invariably 
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introduces undesired decoherence of their quantum states. A master equation model of the quantum evolution of 
one or more qubits interacting with their environment allows one to construct transparent general formulas and is 
therefore very suitable to give both qualitative and quantitative insight into the dynamics of these complex quantum 
systems. 

In this paper we present an analytical solution of a large set of coupled master equations that describe the quantum 
evolution of a particular condensed-matter system, namely the time evolution of an entangled electron spin pair in a 
double quantum dot nanostructure. Even though our model applies to this specific quantum system, the presented 
method of solving the master equations is general and can be applied to study the dynamics of many other complex 
quantum systems. The time evolution of the electron spins is governed by several coherent and incoherent processes, 
each of which depends on time in a simple way as either oscillatory (cosine) or exponential functions. The solution 
we obtain shows how these simple ingredients combine to describe the evolution of the entangled spins in a complex 
nanostructure which consists of several coherent parts. It can be used to predict the occupation probability of all 
quantum states at any given time and to provide analytical estimates of the important time scales in the problem, 
such as the time at which decoherence of the entangled pair becomes substantial^. 

The paper is organized as follows. In Sec. ^] the quantum nanostructure and the assumptions made are described. 
Section IIIII contains the master equations and their solution, with technical details given in Appendices ^ and ^ A 
summary of the results and their range of applicability is presented in Sec. IIVI 

II. THE DOUBLE QUANTUM DOT NANOSTRUCTURE 

The system we consider consists of a double quantum dot nanostructure, which is occupied by two entangled 
electron spins and operated as a turnstile. We studied this system in an earlier paper as a suitable set-up for the 
detection of entanglement between electron spinal. Here we focus on the dynamic evolution of the electron pair in 
the system, which is depicted in Fig. ^ 



(a) (b) 
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FIG. 1: Schematic top view of the double quantum dot nanostructure as discussed in Sec. UTI 
In detail, the structure consists of two adjacent quantum dots in a parallel magnetic field BzZ which are connected 
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to two quantum point contacts (QPCs) via empty quantum channels. A quantum dot is a small metallic or semi- 
conducting island, confined by gates and connected to electron reservoirs (leads) through quantum point contacts. 
If the gates are nearly-closed and form tunnel barriers, the dot is occupied by a finite and controllable number of 
electrons which occupy discrete quantum levels, similar to atomic orbitals in atoms^*^. In our system, the gate between 
the two dots is assumed to be initially open and the dots are occupied by two electronsii [Fig. ^a)] in their lowest 
energy state, the singlet state^^. The gate between the two dots is then adiabatically closed, so that the electrons 
become separated and one dot is occupied by an electron with spin-up and the other by one with spin-down. The 
two spins do not interact anymore and are independently rotated by electron spin resonance (ESR) fields [Fig.^b)]. 
The latter are oscillating magnetic fields which, if the frequency of oscillation matches the energy difference between 
the two spin-split single-electron energy levels, cause coherent rotations of a spin between these levels, analogous to 
Rabi oscillations in a two-level atom. After spin rotation, the electrons are emitted into empty quantum channels by 
opening gates L and R [Fig.^c)] and scattered at quantum point contacts QPC 1 and QPC 2. In a parallel magnetic 
field and for conductances Gqpci(qpc2) ^ /h these QPCs are spin-selective^-, transmitting electrons with spin-up 
and reflecting those with spin-down [Fig. ^d)]. The transmitted and reflected electrons are separately detected in 
the four exits. 

In the next section we analyze the dynamics of the two spins from the moment they are separated and each occupies 
one of the two dots, until both have been detected in one of the four exits. We use a master equation approach in 
which the effects of relaxation and decoherence are included as phenomenological decay ratesi. The solution presented 
is exact under three assumptions: 

• The time evolution during ESR in the dots is decoupled from the time evolution in the channels and exits. 
Physically, this means that the gates between the dots and channels are closed during the ESR rotations, so no 
tunneling occurs out of the dots during that time. 

• Once the electrons are in a channel they cannot tunnel back into the dots, i.e. backreflection of the electrons to 
the dots during their journey to the detectors is neglected. This corresponds to ballistic transport through the 
channels. 

• Once the electrons are in one of the exits they cannot return to the channels, i.e. the electrons are immediately 
detected and absorbed into the detectors. 

III. THE MASTER EQUATIONS AND THEIR SOLUTION 

In the set-up as depicted in Fig. ^ each electron is assumed to be either in a dot, in a channel or detected. This 
leads to a set of 36 possible quantum states represented by a 36x36 density matrix p{t). This set consists of all 
possible combinations AcBa' , with Ag{Z?,C, X} and (7G{t, j} indicating rcsp. the position (Z?=dot, C=channel 
and A"=exit) and the spin direction along z of the electron which started out in the left dot, and Be{D, C,X} and 



(T'G{t, i} representing the position and spin direction of the electron which started out in the right dot. The set is 
given by: 

{i^T^T, ^T^i, DiDt DlDi, CTi^T, C^Dl, CID^, CiDi, D^C^ 
D^C'i, DiC^, Did, CTCt, CTCi, CiCT, CiCi, X^D^, X^Dl 
XiDtXiDi,D^Xt D^Xi, DiXt DiXi, X^Ct X^d, X [Ct 
xia, c^xt c^xi, axt CiXi, X^Xt X^Xi, XiXt XiXi}. (1) 

We number the states in set ^ by the numbers 1 to 36, so l=D'fD'f, 2=D'fDJ, etc. The states labeled by C and X 
do not refer to individual quantum states in the channels and detectors, since in a channel many longitudinal modes 
exist and the detectors consist of many quantum states which form together a macroscopic state. What is meant by 
the states C and X is the set of all channel modes resp. of all quantum states in the detectors. These states thus 
describe the probability of an electron to occupy any one of these channel modes or detector states. We come back 
to why this definition is useful and appropriate in the paragraph below Eq. 0. For long times, the only states that 
are occupied are 33-36, in which both electrons have entered into an exit and the channels and dots are empty. 
The time evolution of the density matrix elements pnm (t) is given by the master equationsi: 

Pn{t) = -^[n{t),p{t)]^^+J2iWnmPm{t)-W^r^Pn{t)), (2a) 
1 

Pn,ra{t) = " , ,„„ " Km Pn,m ) 11^171 (2b) 

for n,m {1, . . . , 36}. The Hamiltonian Ti.{t) describes the coherent evolution of the spins in the quantum dots due 
to the ESR fields and is given by, for two oscillating magnetic fields B^l cos{ujt) x and B^r cos{ujt) x applied to the 
left and right dots respectively, 

T-L{t) = Ho - ^9*Pb cos{ujt) ^ (B^m + ^Bxn)^^m- (3) 

M.Ne{L,R} 

Here TYq is a diagonal matrix containing the energies £"„ {n — 1, . • . , 36) of each state, g* the electron g- factor, /i^ 
the Bohr magneton and &xL(R) ^ 36x36 matrix with elements {(JxL{R))ij — 1 for each pair of states that are 

coupled by the oscillating field BxL(r) and zero otherwise. For g* < the 4x4 upper left corner Hdots{t) of Hit) is 
then given explicitly as 

^ El hARLCos{Ljt) hALRCos{ujt) ^ 



'^dots(0 



HArl COs{ujt) E2 hALRCOs{ujt) 

hALRCos{ujt) E3 TiArl cos{ujt) 

y hALRCos{ujt) fiARLCos{u}t) E4 J 

with El = 2E^ + Ec^ E2 ^ E3 = E^ + Ei + Ec and i?4 — 2Ei + Ec in terms of the single-particle energies E^- and 
El and the charging energy Ec = e^/C, where C is the total capacitance of the quantum dot (assumed to be equal 
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for both dots), Arl ^ An + SL and A^^ = + Sr with Ar(l) = '^''^"af"'^' and 5r^l) = ''''''^2;f"^'^' ■ The 
parameter e, with < e < 1, represents the relative reduction of the field which is applied to one dot at the position of 
the spin in the other dot^. The remaining 32x32 part of the matrix Ti.{t) is diagonal and equal to Ho, since the ESR 
fields are applied when both electrons are located in a dot and the quantum channels do not contain any electrons 
whose spin might otherwise also be rotated by these fields. 

Turning to the transition rates Wnm (from state m to n) in Eqs. (|2a|) . we distinguish between two kinds of transitions: 
1) spin-flip transitions between two quantum states that differ by the direction of one spin only and 2) tunneling 
(without spin-flip) between quantum states that involve adjacent parts of the system, i.e. from dot to channel and 
from channel to exit. The latter are externally controlled by opening and closing the gates between the dots and 
channels. The former are modeled by the phenomenological rate I/Ti.q = Wa^i +^^0x1 with a G {D, C} for spin flips 
in a dot or channel. Here the W's depend on the Zeeman energy AEz = \g*\iiBBz and temperature T via detailed 
balance Wau/Waii; = e^^^/''^'^ , so that 

^"Ti(iT) = ^ l^^-(+|Ag./fe,T ' ae{D,C}. (4) 

The spin decoherence rates Vnm in Eqs. Ij2bp for states n and m with n,m £ {1 • • .4}, i.e. the decoherence rate 
between states in which both electrons are located in a quantum dot, is given by: 

V„™ = + i V (VFj„ Wj™) n,me {!,..., 4}, (5) 

where the W's refer to tunnel rates out of a dot. The coherence between state n and m thus not only depends on the 
intrinsic spin decoherence time T2,d which is caused by e.g. spin-orbit or hyperflne interactions in the dots^**, but is 
also reduced by the (incoherent) tunneling processes from dot to channel^'"^. Similarly, Vnm for all other states n and 
TO is given by 

_ / ife + ^ Ej^„,™(W^j« + Wjm) n, TO e {1, . . . , 16}, but not both n, to e {1, . . . , 4} 
I 00 n e {17, . . . , 36} and/or to e {17, . . . , 36}, 

with the W's tunnel rates from a channel to an exit. Note that energy relaxation processes between different modes 
in the channels, i.e. between modes that contribute to the same set of channel states C, do not affect the transition 
rates Wnm and decoherence rates Vnm for the states where either n or to or both refer to a channel state. The reason 
for this is that these rates refer to resp. spin flip and spin decoherence processes, which are not affected by orbital 
(energy) relaxation and decoherence^^. Hence our definition of the channel states as sets of all modes with the same 
spin does not interfere with the definition of spin relaxation and decoherence of the quantum states. 

With the above ingredients, the coupled equations (|2Jl can be solved analytically. We proceed in three steps: ESR 
applied to the left dot, ESR applied to the right dot and the time evolution after the gates to the quantum channels 
have been opened. During each step only part of the quantum states are evolving in time, while the others remain 
unchanged. This simplifies the procedure to obtain an analytical solution. 
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A. Step 1: ESR applied to the left dot 



Initially, at time t — 0, both spins are assumed to be in the singlet state in the quantum dots, so 

P2(0) = P3(0) = l/2; Pj{0) = Vje{l,4,5,...,36} 
P2,3(0) = P3,2(0) - -1/2; p,j{0)^0 otherwise. 



(7a) 
(7b) 



During ESR applied to the left dot quantum states p^it) — pseit) remain unchanged, since the gates between the dots 
and channels are closed. The coherent evolution of pi{t)-p4^{t) is then governed by the Hamiltonian 

^ El hdLCos{uit) hALCOs{ujt) ^ 

hS^ cos(aji) E2 fi^L cos(ajt) 

TiAl cos(ijjt) E2 MLCOs{u!t) 

y HAl cos{iLjt) hSLCOs{ujt) E4 I 

Including spin-flip rates W^iix and Wdi] and the decoherence rate F = 1/T2^d for both dots^^ we then obtain from 
Eqs. (PI the master equations 



'HESR{i) 



Pi 
h 
h 

Im/5i,2 

Imp24 
Imp34 
Repi^ 
Rep2^3 
Repi^2 
Repi,3 

Rep2,4 
Rep3,4 
Im/5i,4 
Imp2,3 



-8iXveLp\.i - AjJveLpx.z - '2Wdi]Pi + Wd]i{p2 + Ps) 

dLlmpia - ALlmp2,4 + Won + (W^r>iT ~ ^diOpi ~ (2W^DTi + Wdi'\)P2 - WduP3 
-(5Llm/53,4 + ALlm/5i,3 + Won + i^Di^ ~ W^DTi)Pi ^ W^r>TiP2 - {2WD\i + '^Di\)p-i 
-^(P2 - Pi) + (Repi,4 - Rep2^3) - rimpi,2 
-^(P3 - Pi) + Y (Repi,4 - Rep2,3) - rimpi,3 
--^{p^ - P2) - Y (R''3pi,4 - Rep2^3) - rimp2,4 

- Y(P4 - Ps) - y^ (Repi,4 - Rep2^3) - rimp3,4 

- Y (Inipi^3 - Inip2,4) - ^ (Impi,2 - Imp3,4) - rRepi,4 

Y (Inipi,3 - Imp2,4) + ^ (Inipi^2 - Im/53,4) - rRep2,3 

_Ai 
2 

2 

2 

Ai 
2 
<5l 



(Impi4 + Imp2^3) - rRe(5i^2 
(Impi.4 - Imp2,3) - rRepi^3 
(Impi 4 - Imp2,3) - rRep2,4 
(Impi,4 + Imp2,3) - rRep3,4 



(Repi,3 - Rep2 



A, 



(Repi,2 - Rep3,4) - Tlmpi, 



<5l 



(Repi,3 - Rep2,4) 



Al 



(Repi,2 - Re/03,4) - rim/5; 



■2,3j 



(8b) 
(8c) 
(8d) 
(8e) 
(8f) 

(8g) 
(8h) 
(8i) 

(8j) 
(8k) 
(81) 
(8m) 
(8n) 
(80) 



with p,,,(i) = p,j(0 e-*'^* for (ij) e {(12), (13), (24), (34)}, pi,4(i) = Pi,4(i) e^^-* and p2,3(^) = P2,3(0- Eqs. © are 
valid on resonance, so huj = E2 — Ei — E4 — E2 — AEz and within the rotating wave approximation (RWA)^^. p4{t) 



is given by p4{t) = 1 - - p2(0 - Psit)- 

Eqs. (|HJ) can be split into two sets of coupled equations: Eqs. (I8a|l - (PH) and Eqs. ||5j|)- (|8o|) . The solution of the second 
set is straightforwardly obtained and given by 

Repi,2{t) = \ (Re[pi^2(0) - P3a{0)] cos(ALi) - lm[piA{0) + ^2,3(0)] sin(ALt)+ 

Re[/5i,2(0) + p3,4(0)]) e-" (9a) 

Re[/5i,3(0) + p2,4(0)]) e-" (9b) 

Re[/5i,3(0) + p2,4(0)]) e-^* (9c) 
Rep3,4(i) = ^ (-Re[/5i^2(0) - P3,4(0)] cos(AiO + Im[pi,4(0) + ^2,3(0)] sin(ALO + 

Re[/5i,2(0)+P3,4(0)]) e-" (9d) 
Inipi,4(t) = i (Im[pi,4(0) + P2,3(0)] cos(Ait) + Im[pi^4(0) - ^2.3(0)] cos(^it)+ 

Re[/5i,2(0) - P3,4(0)] sin(Ait) + Re[pi,3(0) - p2,4(0)] sin(,5Lt)) e"" (9e) 
Inip2,3(0 = I (Im[pi,4(0) + P2,3(0)] cos(Ai<) - Im[pi^4(0) - ^2,3(0)] cos(5ii)+ 

Re[pi,2(0) - P3,4(0)] sin(ALt) - Re[pi,3(0) - p2,4(0)] sin(<5Lt)) e"". (9f) 



Repi,3(t) = ^ (Re[pi,3(0) - /52,4(0)] cos((5Lt) - Ini[pi,4(0) - ^2,3(0)] sin((5Li) 



Rep2,4(t) = ^ (-Re[pi,3(0) - p2,4(0)] cos(<5Li) + Im[pi,4(0) - p2,3(0)] sm{5Lt) 



In order to solve the set of equations H8a|) - (|5i|l we express pi-ps, Impi_2, Inipi.3, Inip2.4, Imp3^4, Repi^ and Rep2,3 
in terms of new variables xi-xg as follows: 







X2{t)~ 




P2(.t) 




X2(t) + 




P3{t) 


= li-Mt) 


+ X2(t)- 


+ X3it))e- 


Impi,2(i) 


= l{xi{t) + 


XQ{t))e 


-rt 


Impi,3(i) 


= l{x5{t) + 


X7{t))e 


-rt 


Imp2,4(i) 




xr{t))e 


-rt 


ImP3,4(t) 




X(i(t))e 


-rt 


Repi,4(0 


= a;8(i)e"" 






Rep2,3(i) 


= {-xs{t) + 







with Z = Re[pi,4(0) + P2,3(0)]. The transformation pUI) originates from pairwise adding and subtracting those 
equations among 18d|l - l|8i|) which share a common term on the right-hand side, e.g. the equations for Impi.2 and 
Imp3^4. The definition of xi-xg then naturally arises. Physically, the new variables xi, X2 and X3 can be interpreted 
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as Xi(^2) = the probability for the spin in the left (right) dot to be up, and xt, = the probability for the two spins to 
be anti-parallel, each modulated by the exponential dependence on the decoherence rate F. Using (|1U|I . Eqs. (|8a|l - (j5H) 
are rewritten in terms of Xi{t)-xg,{t), which leads to 3 sets of coupled equations. These equations and their solution 
are given in Appendix Eqs. (|10|l at time t — ti, where ti is the time during which the ESR field is switched on, 
thus represent the density matrix elements for the double dot states after the ESR rotation applied to the left dot. 



B. Step 2: ESR applied to the right dot 

Eqs. can also directly be used to obtain the solution after the second ESR rotation applied to the right dot, 
by substituting — > 5jf and Sl — > ^b. in Eqs. © and (jA2|) . and by exchanging xq ^ xj in Eqs. ljA4|l . using pi{ti) 
instead of pi{0) etc. as initial conditions. In order to illustrate this solution, let us consider the initial condition of 
a singlet in the double dot [Eq. (j?))] and let t2 be the duration of the second ESR rotation. In case of no dissipation 
(all VF's — 0) and no influence of ESR applied to one dot on the spin in the other dot we then obtain from Eqs. H1U|I . 
[2|) and ljA4|) for e.g. the occupation probability ^2(^1 +^2) the expression: 



P2itl+t2) = i 



1 + ■{ [cos(OAii) + — s— sin(riAii)] [cosiQAh] + —r— sin(riAi2)] + 
^sin(A-ti)sin(fiAt2)e"**i 1 e~5(*i+*=^) 

riA J 



(11) 



with CIa = i\/4A2 — and A = A^ = Afj. In the absence of decoherence (F ~ 0) the expressions for p2iti + ^2) 
and the other density matrix elements simplify to: 

Pi{ti+t2) = /94(ti + t2) = ^(1 — cos^i cos6'2 — sin(?i sin02) (12a) 

P2{ti+t2) = /93(ti +t2) = ^(1 + cos6liCos6'2 +sin6'isin6l2) (12b) 

Pi,2{ti+t2) = /92,4(ii + ^2) = -^(cos6'i sin6'2 - sin6'i cos6'2) (12c) 

Pi,3(ii+<2) = P3,4(ti + <2) = ^(cos6'i sin6l2 - sin6'i cos6'2) (12d) 

Pi,4{ti+t2) — — -(1 — cos(?i COS02 ~ sin^i sin02) (12e) 

^2,3(^1+^2) = -i(l + cose'icos02 + sin6lisin6l2), (12f) 

with 9i = ilAti and 62 = ^At2- Eq. Hll|l is plotted in Fig. |21as a function of the amount of decoherence F. 

Already for moderate amounts of decoherence T{ti + 12) = 0.001 the occupation probability has become 0.01% less 
than its value in the absence of decoherence P2^'^{ti + ^2) = 0.481 for the set of parametes chosen in Fig. |21 This 
increases to 0.1% for T{ti + ^2) = 0.01. 
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0.49 




FIG. 2: The occupation probability P2(ii +^2) [Eq. lllll l of the quantum state Df D j as a function of the amount of decoherence 
[in units of + 12)]. For r=0 p2 is given by Eq. Ill2bl . Parameters used are A ti = tt/4, A-t2 — tv/8 and all W's=0. 

C. Step 3: Time evolution after the gates to the channels have been opened 

We now turn to the next step in the evolution of the entangled pair in Fig. Q namely the time evolution of the 
density matrix elements after the ESR rotations are completed and the gates to the quantum channels are opened, 
see Fig.^c). From this moment onwards the coherent evolution due to the first term on the RHS of Eqs. ^ stops 
and the time evolution of the matrix elements is solely determined by decay and decoherence rates represented by the 
second terms on the RHS of Eqs. The off-diagonal elements Pi,j{t) then rotate with {Ei — Ej)/h and decay with 
rate Vij: 

PrAt) = P^AtESB.) e'(^»-^^.)(*-*---)/'' ^~V.,it-t.sn) fo^ t > IeSR, (13) 

where Iesb. = ti +^2 and Vij is given by Eq. for i, j G {1 . . .4} and Eq. JBJ otherwise. The initial values Pi.j{tESR) 
for i,j e {!,..., 16} are given by Pm,n{tESR) for m, n £ {1 ... 4} [Eqs. (jTU)) ] with the correspondence in indices 

e {l,5,9,13}^m(n) = l 
e {2,6, 10,14}^m(n) = 2 
e {3,7, ll,15}^m(n) =3 
e {4,8, 12,16} ^ m(n) =4. 

In this way the coherence at time Iesr between any pair of states i,j G {1, . . . , 16} is given by the coherence at Iesr 
between those dot states m,n S {1, . . . , 4} which can (eventually) coherently evolve into i and j, i.e. the dot states m and 
n which have the same spin states as i and j respectively. So, for example, pc^ci.DiciitESR) — PD'\D'\,DiDi{iESR) = 
Pi.iitESR)- Note that Pij{t) — for those states in which at least one electron has reached a detector {i G {17, . . . , 36} 
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and/or j € {17, . . . , 36}), since for those states Vij = oo. This corresponds to the assumption of immediate detection. 

In the remaining part of this paper we focus on the evolution of the populations pi{t)-p3e{t) for times t > Iesr 
under the following conditions: 

• We neglect the possibility of spin flips in the dots, i.e. we set Wdtx — Wdii — 0- This based on the fact 
that Tijj is known to be much longer (0.85 ms at magnetic fields = ST)^^ than the time required to travel 
through the channels to the exits. This assumption is not essential to obtain an analytical solution; it only 
simplifies the resulting equations. 

• We assume that the tunnel rate Wt out of the dots into the channels is equal for spin- up and spin-down electrons, 
i.e. the two electrons tunnel out of the singlet state with a negligible time delay tdeiay in between, and that spin 
is conserved during this tunneling process. TypicalljiSS tdeiay ~ 10^^^ s, which is much less than the travel time 
through a channel ~ 10^^'^s. 

• The tunnel rate We through the QPCs is taken to be constant and equal for spin-up and spin-down electrons, 
i.e. the set-up is assumed to be constructed in such a way that the detection time for spin-up and spin-down 
electrons once they have reached the QPCs is the same. 

• Spin flips in the exits are neglected, i.e. detection is assumed to be faster (with typical times ~ 10"""^^ s)2i than 
the spin-flip rate (3> 10^^^ s)^ in the detectors. 

The evolution equations for pi{t)-p^Q{t) for times t > tESR Eire then given by the master equations 



p, = -2WTP^ for ie{l,...,4} (14) 

Pi = WrPi-i + Wc]iPi+2 - {We + Wt + Wci^)pi i e {5, 6} (15a) 

P, = WTP^-4 + Wci]Pi~2 - {We + Wt + WcM)pi i e {7, 8} (15b) 

P, ^ WTP^-S + WcllPi+l-{WE + WT + Wcn)p^ iG{9,ll} (15c) 

P^ = WTP^-8+Wcl^p^-l-{WE + WT + Wcn)p^ «e{10,12} (15d) 

Pi3 = Wt{p5 + P9) + Wcn{Pi'i + Pi5)~2{We + Wcii)pi3 (16a) 

pi4 = Wt{p6 + Pia) + Wci^pi3 + WcnPie ^ {'^We + Wcn + ^cii:)pi4: (16b) 

Pi5 = Wt{p7 + pii) + Wci^pi3 + WcnPi6 - {"^VVe + Wcu + ^ci'i)pi5 (16c) 

Pw = Wt{p8 + P12) + Wci\{pi4 + P15) -'2{We + Wcu)Pw (16d) 



Pi{t) = WEPi-12{t - ttravel + tESR) " Wt P^{t) i G {17, . . . , 24} 



(17) 
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ft W = Wrft-sW + WEP^-12{t - ttravel + + W^CTift+l W " i^E + Wcl^)p^{t) 

fori €{25,27} (18a) 

ft(t) = WTP^-s{t) + WEP^-12{t " Uravel + tESu) + Wcl1P^-l{t) - {We + Wcn)Pt{t) 

for {26, 28} (18b) 

p,{t) = WTP^-8{t) + WEp^~l&{t " Uravel + tESIi) + WcuP^+2{t) - {We + Wci^)p^{t) 

fori €{29,30} (18c) 

p^{t) = WTP^-&{t) + WEP^^le{t " Uravel + tESR) + Wci]P^-2{t) - {We + Wcu)P^{i) 

for {31,32} (18d) 



P^=WE{p^^B+ P^^A) ^ G {33, . . . , 36} . (19) 

Here ttravel > tESR denotes the earliest time at which an electron has traveled through the channels and reached 
an exit. For times t < ttravel, Pi{t) — P36{t) is thus given by Eqs. (|14|l - p9|) for We — 0, since at those times no 
electron can have arrived at a detector yet. The above sets of coupled equations can be solved one by one: first those 
for pi{t)-p4{t), then once the latter are known those for Pn{t)-pi2{t) (hi the pairs (5,7), (6,8), (9,10) and (11,12)), 
then pi'i{t)-piQ{t) and pn{t)-p2i{t), subsequently P25{t)-p32{t) (in the pairs (25,26), (27,28), (29,31) and (30,32)) and 
finally P33{t)-p3e{t). Proceeding in this order and using initial conditions 

f Eqs. Uni) for i = 1, ... ,4 

P^{tESR)^{ (20) 

[ for i = 5, . . . , 36 

we obtain for pi{t)-p4^{t), the states in which both electrons are located in a dot: 

P^{t) = p^{tESR) e-2M'T(t-tBSH) j £ { 1 , . . . , 4} , i > tESR- (21) 

Next, we find for P5{t)-pi2{t), which correspond to the quantum states in which one electron is located in a dot and 
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the other in a channel, from Eqs. 115(1 : 

Psit) = A5,7,i^3e-^^^^(*-*^^«) +55^7,1,3 e-(^^+'^^)(*-*^^«) +Ci,3e-2'^^(*-*«s«) (22a) 

p^(^t) = -As 7 1 3e->^--c(*-*-s«) + ^^^^5 7 1 3 e-(>^-+^-)(*-*-s«) + 3 e-2W'T(t-*BSH) (22c) 

P9(0 = A9,io,i,2e-'^^^^(*-*^^«^ +5940,1.2 6"^^^+'^^^'*"*^^'^^ +Ci,2e-2'^^(*-*^^«) (22e) 

Pio(t) = -Ag 10 1 2e-^--^(*-*-^«) + ^^^^9 10 1 2 e-(>^-+^^-)(*-*-««) + Di 2 e-2^-(*-*-^«) (22f) 

pii{t) = Aii,i2,3,4e-^^^'=^(*-*^^") +Bii42,3,4e-(^^+'^^)(*-*^^") +C3,4e-2'^^(*-*^««) (22g) 

Pi2(t) - -All i2 3 4e-'^^^^(*-*^^«) + ^^^Bii i2 3 4 e-('^^+^^)(*-*^^«) +i:'3 4e-2'^^(*-*^^«', (22h) 

Wen 



where 



Wetc = We + Wt + Wen + Wen (23a) 

^ ^ W^CiTP»(^ggfl) - WcnPlitESR) WT{-WcnPk{tESR) + W^CTiPi(^i;Sfl)) /2ou^ 

'''''''' " W^CTi+W^CiT (W^CTi + W^CiT)(W^s-W^T + M^cu+W^CiT) 

-Bij,fc,i = ^^^^ _^^^^^ [AftgSfl) + Pj{tESR) - _ ^^ {PkitESn) + Pl{tESR))] (23c) 

^ Vt^T (VFg - T^T + W^cTi)p»fti;5fl) + WcnPj(tESR) .23dl 

w^ij -Wt We-Wt + Wen + Wen 



^ ^ Wt WcnPii^ESR) + {We - Wt + Wcn)Pj{tESR) . 
"'^ We ~Wt We-Wt + Wen + Wen ' ^ 



For times t < ttravei, the evolution of P5{t)-pi2{t) are given by Eqs. H22|) with We — 0. For times t > Uravei these 
populations are given by Eqs. H22|) with Iesr ttravei- 

In order to obtain the solution for pi^{t)-pi(,[t), which corrresponds to the situation in which both electrons are 
located in a channel, we rewrite the equations for pi^-pi^ as 

Pi3 = WT{pb + Pq) + Wen{pu^' P15) -'iiWE + Wen)Pi^ (24a) 

P14 + P15 = Wt{pq + P7 + pw + Pii) + '2WcnPi3- {'^WE + We]i + Wen){Pu + Piz) + '^WenPi(i (24b) 

P16 = Wt{p9. + P12) + Wen{Pi'i + P15) -"^-iWE + We]i)piQ (24c) 

P14-P15 = Wt{p6 - P7 + Pw - Pii) - {"^We + We]i+Wen){Pii- Piz)- (24d) 

Eqs. 124|) consist of 3 coupled equations H24a|l - H24c|l and a separate one, Eq. (|24d|) . We first solve the latter and then 
the first three. In each case the solution is a combination of a homogeneous and a particular solution. Taking from 
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i 


A, 






17 


-W_bA5_7_i 3 


-65,7,1,3 


We r< 


18 


— W£;A6 8,2,4 


-06,8,2,4 




19 


WeAb,7,1,3 
We+2VKc 


55,7,1,3 


3 


20 


^6,8,2,4 

TIT" 1 '^1^/' 

W E+-^vV c 


-06,8,2,4 




21 


— WeA9,10,1,2 
We+2VKc 


^9,10,1,2 




22 


^-5^9,10,1,2 
We+2Wc 


^9, 10, 1,2 


M/y -^1.2 


23 


— W^E ^11,12,3,4 


Bll,12,3,4 


We 


We+2Wc 


24 


^11,12,3,4 
We+2Wc 


^11,12,3,4 


We n 



TABLE 1: Coefficients Ai, and d in Eqs. TO 

now onwards W^cTi.=^CiT = W^c^^ we obtain: 
Pi3(i) - 



^g-2(H'E+Wc)(t-tBSH) I l^g-2WE(*-tESH) _ 1^ g- 2 (IVe +2Wc ) (t-tESfl) 

2 2 



Pi4(i) 



+ His e"(^^ + ^^+^^^)^*"*^^^^ + Kl3 e-i^E + WT){t~tESR) ^ ^^^^ g-2WT(t-tESfl) 

_ pj „-2{WE + Wc){t-tESR) 1 1 p p-2WE(t-iBSJi) I }_p „-2{WE+2Wc){t-tESR.) 

2 ' 2 

+ i7l4 e"(^^ + '^^+2^'=^)(*~*^Si?) _|_ g-(WE + WT)(t-tESi?.) _^ g-2Wr(t-tBSB) 

^ _ 771 -2(VKE+lVc)(t-tESH) I }_p -2WE{t-tESR) 1 1 p -2(H'E+2H'c)(t-tESH) 

2 2 

^ ^~2{WE + Wc){t-tESR) I l^g-2WE(t-iESH) „ 1 ^ g-2 (IVe +2Wc ) (t-tESfl) 

2 2 



(25a) 



(25b) 



(25c) 



(25d) 



The coefficients in Eqs. (I25|l are given in Appendix IbI Also here, Pi3{t)-pie{t) for times t < ttravei are given by 
Eqs. (|25|l with We — 0, and for times t > ttravei these populations are given by Eqs. (|25|l with tssR ~^ ttravei- 

The solution of the next set, pn{t) — P24(i), corresponding to the states in which one electron is located in a dot 
while the other has reached a detector, is given by 



Pi{t) ^ Ai e~(^^ + ^^ + 2M'c)(t-tt™„.!) _ g-(W'E + WT)(t-ttra,,el) _ (J. g " 2 (t " * t,a„.i ) 

[p^{ttravel) ' A, + B, + C^] e-M'^(*"*— ') for i g {17, . . . , 24}, t > ttravei 



(26) 



and pi{t) = for t < ttravei- The coefficients Ai, Bi and C,; in Eqs. H26|l are given in Table |I| Next, we solve 
for P25{t)-p32{t), the states in which one spin has reached a detector, while the other is still in a channel, in the 
pairs pi{t)&zpj{t) £ {p2b{t)^P2dt), P27{t)kp2s,it), P29{t)kpzi{t), and pso (0-^/032 (i)}, see Eqs. (UHl. For each pair the 
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solution is given by, for times t > ttravei- 

M,,5e-(^-+'^-)(*-*— ) +Af,,6e-2^-(*-*— +M,,7e-^-(*-*— ) (27a) 

M,,5 e-(^-+'^-)(*-*— ) + M,, 6 e-2^-(*-*— ) + M,,7 e-^^^*"*— \ (27b) 

and pi{t) — pj{t) = for t < ttravei- The coefficients Pij, Qij and Mi^i . . . Mij for i,j e {25, . . . , 32} are given in 
Appendix IbI 

Finally, we obtain the time evolution of the states pssit) — pseit) in which both electrons have reached an exit. This 
is given by, for times t > Uravei , 

' - "^^X We We + 2Wc 

2(iy£; + iyc) 

^^m,3 + Mn^'i -2{WE+2Wc){t~tt,-a^,i) , -^7n,4 + Af»,4 -(WE+WT+2M'c)(t-tt™„el) 

2(1^^ + 2Wc) We + Wt + 2Wc 

W^B + Wt 2Wt 

Mm J + M„,7 vyT(t~tt.a„.i) 

+ iy_B( sum of all previous coefficients, so "''^ — — + '^'"'^ + ■ • • ) (28) 

We We + 2Wc 



for 

(j, m, q) G {(33, 25, 29, 26, 31), (34, 26, 30, 25, 32), (35, 27, 31, 28, 29), (36, 28, 32, 27, 30)}. 

Special case. In order to illustrate the solution p8() . we now derive explicit expressions for pssit) and p34(t), the 
probabilities that a spin-up is detected in the left detector and resp. a spin-up or a spin-down in the right detector, 
for the special case of F = (no decoherence in the dots) and Wcu = W^ciT — (no relaxation in the channel). This 
corresponds to the situation in which the time evolution occurs in the absence of any decoherence and dissipation 
mechanisms in the dots and channels and only depends on Wt, the tunnel rate from dot to channel, and We, the 
tunnel rate from channel to exit. 

We are interested in finding P33{t) and ps^it) for times t > ttravei [since P33{t)=p34(t)—0 \f t < ttravei]- To that end, 
we first calculate Pj (ttravei) for j < 16 from Eqs. ()21|l . H22|l and H25|l and then all coefficients entering the expressions 
for P33W and P34(i) in Eqs. 
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For pi{ttravei)-pi6ittravei) we then obtain: 

Plittravel) = P4 (^tra^e; ) = Pi (^BSfl) e""'^^ "^^^^-^ (29a) 

P2{ttravel) = P3 (W.ei ) = P2 (iiJSfl) c"'"^^ (29b) 

P^{ttravel) = (t^s^) e^'^^ ' (1 - e"^^ ' ) ^ G {5, 8, 9, 12} (29c) 

P^{ttravel) = p2 (tfisfl) e"^^ ' " e"'^^ "'^^''^ ) {6, 7, 10, 11} (29d) 

Pl3(it™.e;) = Pl6(W.ei) =Pl(^BSi?j(l-e-'^^(*"— '-*^^«))2 

Plittravel) = Pl5 (W.ei) = P2 (^BS/?) (1 " e"*^^ )2 . 



(29e) 
(29f) 



Eqs. ((23) form the initial conditions that appear in the expressions for P33{t)-p3Q{t) [Eqs. (j^HH- We then find V 

i ^ ^travel- 

P,3{t) = [^OuiUravel] - ^^^^^^ P^Uravel) + (^^^^t)^ PlUtra.el)) e-^'^-^*— ) + 

(-2pi3(t*.a.e.) - ^^^^^P5(W.eO + -^^.^^ Pi (t*™.eO) e^^^ ) + 

-(P5(itra«ei) + Pl(W.ei)) e"*^^ f*"*"— ' ' + 



2We 
We-Wt' 



■pi(tt™.eOe-^'^-(*-*— ) + 



Pl3{ttravel) + "2 Psittravel) + Pl{ttravel) (30a) 
P34(t) = (pi4(t..a.eO - ^^'^^^^ P.iUra.el) + (^^^^t)^ P^l^^^^^eO) e-^>^-(*— ) + 

(-2pi,{Ura.el) ^ ^^^^-^^ P6 ( W.e.) + P2 (t*.a.eO) e-^-(*— ) + 

V yy E ~ yVx yyE~ J 

(Peittravel) + PiiUravel)) e"'^^ ^ + 



We - Wt 

2We 
We-Wt' 



■p2{ttra.el)e-'^-^'-'^--^ + 



{We - WtY 

Pli{Uravel) -\- 2 pQ{ttravel) + P2{ttravel)- (30b) 



One can see directly from Eqs. (|30|l that the time dependence of ^33 and P34 is determined by five exponential 
functions, whose relative magnitude depends on the ratio between We and Wt- This is illustrated in Fig.|21 which 
shows Eqs. (|30ll as a function of t — ttravei for various rates We and Wt- For Wt ^ We the time needed to reach 
the stationary state (the average detection time) is dominated by the term ~ e~'^^'-*~**™''='\ whereas for Wt ~ We 
the terms e~^'^^(*"**'-<""='\ e-f'^^+^^^f*-*'™"-') and e-^w'rlt-tt™^.;) dominate. 
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0.4 



0.2- 



0.3- / 



0.1-1/ 



0, 








1 



2 



3 



4 



5 



^ ^ travel ^^^^ 



FIG. 3: The probabilities pss to measure two spin-up electrons and P34 to measure a spin-up and a spin-down electron in 
the left and right exits resp. for times t > Uravei- Parameters used are Oi = n/2, O2 = tt/S [so that pi(fBSfl)=0.154 and 
P2{tESR)=0M6], ttravei - tssR = 0.1 ns, We = 10^° s"^ and Wt = 10''s"^(9.9 lO'' s"^) for the solid (dashed) curves. 



In summary, wc have presented an analytical solution of a set of coupled master equations that describe the time 
evolution of an entangled electron spin pair which can occupy 36 different quantum states in a double quantum 
dot nanostructure. Our method of solving these equations is based on separating the time evolution in three parts, 
namely two coherent rotations of the electron spins in the isolated quantum dots and the subsequent travel of the 
electrons through two quantum channels. As a result of this separation, the total number of master equations is 
split into various closed subsets of coupled equations. Our analytical solution is the first of its kind for a large set 
of coupled master equations and the same method can be used to study and predict the quantum evolution of other 
quantum systems which are described by a large set of quantum states. This type of analysis complements numerical 
approaches to study the dynamic evolution of complex quantum systems and allows to obtain qualitative insight in 
the competition between time scales in these systems. 
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APPENDIX A: SOLUTION OF EQS. (ISat-dSill 



Using the substitution Eqs. H1U|I . Eqs. H8a|l - l(5H) transform into: 



ii = ~{Wdi^ + WDn ~ ^LX5 + WDn(^ 

X5 = Alxi —e 



rt 



X2 = -{Wdi]' + WDn -^)x2 - Slx^ + WDn^ 

r rt 

XA = OlX2 - —e 



rt 



Al 5l 



with Z = Rc[pi_4(0) + P2,3(0)]. In deriving Eqs. (|A1|) we have used that 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
(Ale) 
(Alf) 

(Alg) 
(Alh) 



Pi = 1 - Pi - P2 - P3 

R.ep2,3 = -Repi_4 + Ze^'"*. 

Equations. IjAip consist of three sets of coupled equations, (jAlap - (jAlb|l . HAlc|) - (|Ald|) and HAle|) - (|Alhl) . The sohition 
of the first two sets is given by: 

-{Wdm + WDi^ - r)(xi(0) - Ai) + 2Ai(a;5(0) - ^2)^ 



Xl{t) 



X5(t) 



X2{t) 



Xi{t) 



/2AL(a:i(0) - Ai) + {Won + Wdi^ - ^(2:5(0) - A2) 



rt 



V 217a 

{X5{0)- A2)cosnAt] e-i'-^^n+Won~r)t _^ 

{Wpi^i + Wpi^ - r)(x2(0) - A3) + 2(5l(x4(0) - Aj) 
2ns 

{X2i0)-A3)cosnst] e-5(M'«Ti+M^i>iT-r)t + j^^^rt 

2SL{x2io) - A3) + [Won + Wdi^ - r){xi{o) - a^) 



sin Q/^t + 



sin Q/^t 



sin Qst 



2n>i 



sin flgt 



(a;4(0) -^4)cos175<] e-5(M'«Ti+M'z,iT-r)t + 



rt 



(A2a) 



(A2b) 



(A2c) 



(A2d) 
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with 



^3 



Aj + 2TWdu 

2{^l +T{WDn+Wni^)) 

51 + 2TWdu 
2{5l+T{WDU+WDi^)) 

2{5l+T{WDn+WDi^)y 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 
(A3f) 



So far no approximations have been made, apart from assuming the decoherence rate T to be equal for aU off-diagonal 
terms of the density matrix p [Eqs. JHl]. In order to obtain the solution of the remaining equations ljAle(l - (|Alh|l we 
assume 5l = (no influence of the ESR field on the spin in the right dot) and = VFdxt = ^^^"^ findSi; 

"r (2:3(0) -i)+2ALX7(0) 



X3{t) 

xe{t) 

X7{t) 



2VL 



A 



sinf^A^ + (^3(0) — -) cos r^At 



X6 (0) cos Ait + (2x8 (0) - Z) sin Ait 

-2Al (x3(0) - \) - . ^ ^ ^ ^ , 
smiZA^ + 2:7(0) cos iZAi 



2n; 



62 ■ 



xsit) = -[-2:6(0)sinAit+(22:8(0)-Z)cosAii + Z] 



(A4a) 
(A4b) 
(A4c) 

(A4d) 



with i7/ 



2 _ p2 



APPENDIX B: COEFFICIENTS OF EQS. (HH) AND HTh 

The coefScients in Eqs. (|25f) are given by: 

— PlsitESR) + PieitsSR) + — HiQ + Ki3 — Kiq + L13 — Liq 



E 
E 
F 
F 
His 



H 



14 



K 



14 



2 

PliitESR) — PlsitESn) — Hi4 + Hi5 — Ki4 + if 15 — L14 + L15 

2 

PlsjtESR.) + PltjtESR) + PlbjtESR) + PwjtESR) " 2(ifi4 + X15) — L13 ~ L14 — L15 - Lie 

2 

— PiajtESR) + Pl4(^£gii:) + Plbi^ESR) " PieitESR) " 2(J?i4 + ifis) + Lia — £14 — Li5 + Lig 

2 

WT[(Tyg - VFt + W^c)(A.7,l,3 + ^,10,1,2) + Wc{Ag^8,2A + ^n, 12,3,4)] 

{We - Wt){We -Wt + 2Wc) 
Wt[VFc (^5,7,1,3 - ^11,12,3,4) + {We -Wt + Wc){Aq^s.2A - ^9,10,1,2)] 



{We ~ Wt){We ~Wt + 2Wc) 
Hi5 = -//i4[(5, 7, 1,3)^ (6, 8, 2, 4), (9, 10, 1,2)^ (11,12,3,4)] 

Hw = -i/i3[(5, 7, 1,3)^ (6, 8, 2, 4), (9, 10, 1,2)^ (11,12,3,4)] 

Ku = HisiA^B) 

Wt[W^c (^5,7,1.3 + ^1142,3,4) + (We -Wt + Wc){Bg^s,2a + 59,10,1,2) 



{We - Wt){We -Wt + 2Wc) 
Ki5 = ifi4[(5,7,l,3)^(6,8,2,4),(9,10,l,2)^(ll,12,3,4)] 

Kie = Hie{A^-B) 

Ll3 = WT[{2{WE-WT + Wc f -W^){Ci^2 + Ci^3) + WS{D2,i + D3,4) + 

Wc{We -Wt + Wc){C2a + CsA + ^1.2 + ^1,3)]/ 
[A{We - Wt){We -Wt + 2Wc){We -Wt + Wc)] 

Ll4 = Wt[Wc{We -Wt + Wc){Ci^2 + Cl,3 + D2A + ^3,4) + 2{We -Wt + Wc?{C2,i + i?l,2) 

+W^(C3,4 - C2,4 - £"1,2 + Di.:i)]/[^{WE - Wt){We -Wt + 2Wc){We -Wt + Wc)] 
Liz = -^i4(C2,4 ^ C3,4 , £'1,2 ^ -Di,3) 

L16 = Li3(Ci,2 ^ -02,4 , Ci, 3 ^ 153,4). 



20 



For (i,j)=(25,26) the coefRcients in Eqs. (|27|l are given by 

7 

2 



^25, 26 = ^26,25 — ~ ^(^^25,fc + -M26,fc 
k=l 

(925,26 = ^Q26,25 = ^(-^^^ 



25,1 



25,4 



25,5 



fc=l 

7 

25, fc ^ ^^26,fej 



2 

fe=i 



(T4^£; + M^c)-B + WcE 



We + 2Wc 

1 

M25,2 = M26,2 =-2^ 

Wv 

M25,3 = -M26,: 



2{We + 2Wc) 

We[-{Wt + Wc)iiWE + 2Wc)Hi3 - WtA5,7,1,3) + WcHWe + 2Wc)Hu ~ WtA,8,2,4)] 



Wt{We + 2V^c)(Wt + 2Wc) 

M26,4 = M25,4(i?13 ^ i^l4,A,7,l,3 ^ ^6,8,2,4) (B2) 
{Wc - WT)iWEKi3 - Wt55,7,1,3) + Wc(WEi^l4 " T^tS6,8,2,4) 



Wt{Wt - 2Wc) 

^j. ^ We[{We - 2Wt + Wc){Li3 - Ci^3) + WcjLu - ^2,4)] 
^^''^ {We-2Wt + 2Wc){We-2Wt) 

-^26,6 = Af25,6(-^13 ^ -^14, Cl,3 ^ C2,4) 

M25,7 = Wt [{We ~Wt + Wc)[-Aij + + C17] + W^c[-^i8 + ^18 + Cig]] / 
[{We -Wt + 2Wc){We ~ Wt)] 

M26.7 = M25,7(17^ 18). 



The coefficients in Eqs. ((23 for (i,j)=(27,28), (29,31) and (30,32) are obtained from Eqs. (HU by replacing indices 
as given in Table Hll 
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(27,28) 


(29,31) 


(30,32) 


25 ^ 28 


25 29 


25 32 


26 27 


26 ^ 31 


26 ^ 30 


13 16 




13 ^ 16 


14 15 


14 15 




17 ^ 20 


17 21 


17 ^ 24 


18 ^ 19 


18 ^ 23 


18 22 


E ^ -E 




E ^ -E 


E ^ -E 


E ^ -E 




^5,7,1,3 ^ —^6,8,2,4 


(5,7,1,3) ^ (9,10,1,2) 


j45,7,i,3 — v4ll,l2,3,4 


^5,7,1,3 56,8,2,4 




-65,7,1,3 — > -Bll,12,3,4 




(6,8,2,4) ^ (11,12,3,4) 


j46,8,2,4 — > —^9,10,1,2 






^6,8,2,4 — > 59,10,1,2 


Cl,3 -D2,4 


Cl,3 ^ Cl,2 


Cl,3 — » -D3,4 


C2,4 Di^3 


C2,4 — * C3,4 


C2,4 Di^2- 



TABLE II: Required substitution of indices and coefficients in Eqs. 1271 in order to obtain tlie corresponding coefficients for 
p,{t)&,p,{t) witli (ij) G {(27, 28), (29, 31), (30, 32)} 
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See the calculation of WcTi VKcxt in Refi-. 

For the general case the expressions are rather lengthy, but not qualitatively different. 

The solutions for xs{t) and xs{t) given here are valid for 5l ~ 0, i.e. assuming the W's to be zero is not necessary. 



